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1. Introduction 

Perturbative quantum field theory is progressing fast. New high-precision experi- 
ments require calculation of higher radiative corrections — multiloop Feynman di- 
agrams. Recently, some calculations have been done which would seem impossible 
only a few years ago. This is due to the high degree of automation of the process 
of generation, analyses and calculation of Feynman diagrams, which is achieved via 
extensive use of computer algebra (see |IJ for review and references) . All such calcula- 
tions are performed in the framework of dimensional regularization 0, i.e. diagrams 
are calculated as analytical functions of the space-time dimension d = 4 — 2e. 

The integration-by-parts method pj was invented for calculation of three-loop 
massless propagator diagrams. It is the most systematic method of those currently 
used, and the most appropriate for computer-algebra implementation. It was first im- 
plemented clS db SCHOONSCHIP £| package MINCER f|, and later re-implemented [g 
in FORM (in fact, FORM was created mainly to run MINCER). Since then, MIN- 
CER has been the engine behind most of spectacular successes of perturbative field 
theory. Some of these calculations, with gigabyte-size intermediate expressions, are 
among the largest computer-algebra calculations ever undertaken. 

Integration-by-parts was used for other classes of problems, too. Many inter- 
esting physical results have been obtained with packages for calculating two-loop 
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on-shell massive diagrams and three-loop vacuum diagrams with a single mass (see, 
e.g., H H [L0||). Reduction of two- loop propagator integrals with generic masses and 



momentum to a finite set of bases integrals has been achieved [11]. First three- loop 



on-shell calculations have been done recently [12, [13]]. 

Several years ago, an interesting new approach to heavy-quark problems in 
Quantum Chromodynamics has been formulated — Heavy Quark Effective Theory 
(HQET), see, e.g., [13, |T3| for review and references. In collaboration with David 
Broadhurst, I applied the integration-by-parts method for calculating two- loop prop- 
agator diagrams in HQET ||16|| . Since then, the algorithm suggested was used in a 



large number of physics applications. A short review of the integration-by-parts 



method as applied in heavy quark physics is presented in |L7]. 

In the present work, I apply this method for calculating three-loop propagator di- 
agrams in HQET. Three-loop anomalous dimensions and spectral densities in HQET 
are necessary for a number of physics applications, such as improved extraction of 
the B meson decay constant from lattice simulations and from QCD sum rules. 
HQET lagrangian does not involve mass in the leading order, and, therefore, the 
problem is quite similar to the massless one. My aim is to produce a reliable pack- 
age for three-loop HQET calculations, which could play the same role as MINCER 
in massless theories. I call it Grinder. Some complication comes from the fact that 
there are two kinds of lines now — massless propagators and infinitely-heavy ones, 
and hence the number of diagram topologies is substantially larger. 

In section |2], we consider three-loop HQET propagator diagrams with one- or 
two-loop massless or HQET propagator subdiagrams. To this end, we first recall 



well-known results for massless B and HQET \TE\ one- and two-loop diagrams. After 



that, diagrams with two-loop insertions, or with two one- loop insertions, are easily 
calculated. Two-loop HQET diagrams with a single one-loop insertion are dealt with 
in a manner similar to the plain two-loop diagrams. In section [|, we consider proper 
three-loop HQET propagator diagrams. Some details of implementation and testing 
are presented in section |j. 

2. Diagrams with lower- loop propagator insertions 

2.1 Two-loop massless propagator diagrams 

The one-loop massless propagator integral (figure P can 
be easily calculated by using the Feynman parameterization or 
Fourier transform to the coordinate space and back: 

u \ u 2 Figure 1: One- 

D\ = —k 2 , D 2 = —(k + p) 2 , (2-1) loop massless prop- 

T(ni + n 2 — d/2)T(d/2 — rii)T(d/2 — n 2 ) agator diagram. 




G{n u n 2 ) 



T{n 1 )T(n 2 )T(d - ni - n 2 ) 
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Figure 2: Two-loop massless propagator diagram 



The integral with a numerator can be written as a finite sum 

T d/2r_ 2\d/2-m-n 2 



P n {k)d d k 



m \ f* / 



X 

J71 



, (2.2) 

k— >p 



. r(ni + n 2 — m — d/2)T(d/2 — n\ + n — m)T(d/2 — n 2 + m) 

G{ni,n2]n,m) = — — 

1 (niji [n 2 )i [d - ni - n 2 + n) 

where P n (k) is an arbitrary homogeneous polynomial: P n (Xk) = X n P n (k). 
We write the two-loop propagator integral (figure |2]a) as 

d d k\d d k 2 j, 2\d-Vnr(/ ^ 

n» 2 n »3 n »4 n » 5 = ("P ) ^ n *G(ni, n 2 , n 3 , n 4 , n 6 ) , 



D x = -k\, D 2 = -k 2 2 , D 3 = -(k 1 +p) 2 1 D 4 = -(k 2 +p) 2 , 
D 5 =-{k x -k 2 ) 2 . (2.3) 

It is symmetric with respect to 1 <->• 2, 3 <-> 4, and also 1 <-> 3, 2 <->• 4. If one of the 
indices is zero, it can be easily calculated using (|2.1| ) (figure 06, c) 

Gf(ni, n 2 , n 3 , n 4 , 0) = G(ni,n 3 )G(n 2 , n 4 ) , (2.4) 
G(0, n 2 , n 3 , n 4 , n 5 ) = G(n 3 , n 5 )G (n 2 , n 4 + n 3 + n 5 - d/2) (2.5) 

(and symmetric relations). 

Applying the operators d\ ■ {k\ — k 2 ) and d\ ■ k\ (where di = d/dki) to the 
integrand of (|2.3|) , we obtain the recurrence relations for G(rii,n 2 , n 3 , n 4 , 715) (known 
as triangle relations ||) 

[d - rn - n 3 - 2n 5 + ni l + (2- - 5") + n 3 3 + (4~ - 5")] G = , (2.6) 
[d-n 3 -n 5 - 2n x + n 3 3 + (l - 1") + n&+(X - 1")] G = , (2.7) 

where, for example, 

l ± G(ni, n 2 , n 3 , n 4 , n 5 ) = G(m ± 1, n 2 , n 3 , n 4 , n 5 ) . (2.8) 
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Figure 3: Basis two-loop massless propagator integrals 

Of course, more relations are obtained by symmetry. Another interesting relation is 
derived by applying the operator ■ (k 2 + p). Substituting the general form of the 
relevant vector integral, we arrive at 

+ n 4 - ni - n 2 - n 5 + (^d - n\ - n 2 - n 3 - n 4 - n 5 ^j (4~ - 2~) + (2.9) 

+n 3 3 + (4~ -5") G = 0. 

This formula was derived long ago by S.A. Larin in his M.Sc. thesis |18[] (again, 
similar relations follow by symmetries). 

If indices of two adjacent lines are non-positive integers, the integral contains a 
no-scale vacuum subdiagram and hence vanishes. The cases with zero indices are 
given by fl2.4Q, (|2.5|). When n 5 < and n 3 7^ 1, n 5 can be raised by (|2.9|); the cases 
r?-! ^ 1, ri2 ^ 1, 714 1 are symmetric. The case n 5 < 0, n\ = n\ = n 3 = n 4 = 1 is 
handled by 

[{d - 2n 5 - 4)5+ + 2{d - n 5 - 3)] G(l, 1, 1, 1, n B ) = 

= 21+(3- - 2-5+)G(l, 1, 1, 1, n 5 ) , (2.10) 



which follows from (|2.6|) and ( p.7|) at rii = ri2 = ^3 = 714 = 1 (note that the terms 
in the right-hand side of (|2.10|) are trivial for any n^; for n*, < 0, they vanish). 
When n 2 < 0, it can be raised by (|2.9|); the cases n\ < 0, n 3 < 0, n 4 < are 
symmetric. 

We are left with the most important situation when all the indices are positive. 
Applying (|2.6| ), we reduce n 2 , n^, until one of them vanishes. Then (|2.4|) and (|2.5|) 
apply. If max(ni,n 3 ) < max(n 2 ,n4), it is more efficient to lower n\, n 3 , 71,5. 

All one-loop integrals (figure |I]) with integer n\^ 2 are proportional to G\ = 
G(l,l), the coefficient being a rational function of d. All two-loop integrals with 
integer indices reduce to G\ (figure ||a) and G 2 = (7(0,1,1,0,1) (figure |3|fe), with 
rational coefficients. Here 

r = 1 r(l + ne)r" +1 (l-6) 

" (n+l-n|) B ((n + l)|-2n-l) B r(l-(n+l)c) ' l ' J 

where (x) n = T(x + n)/T(x) is the Pochhammer symbol. 
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2.2 Two-loop HQET propagator diagrams 

The one-loop HQET propagator integral (figure |) can be 
easily calculated by using the modified Feynman parameteriza- 
tion (see, e.g., [16|, [T^J), or Fourier transform to the coordinate " ~ 
space and back: Figure 4: One- 

f ddk ■ d/2i o \d-2n 2U \ loop HQET propa- 

J WW = m ( } (ni ' n2) ' g ator dia g ram - 

I{n u n 2 ) = — 2.12 

r(ni)r(n 2 ) 

(here v is 4- velocity of the heavy quark, v 2 — 1, and co> = p ■ f is the residual energy). 
Similarly to (|2.2|) , we obtain 

p n{k)d d k _ d/2/_ 2 \d-2na 



^ /(ni ' n2;n ' m)L ^(-il;^) P " (A;) 



r(ni + 2n 2 - n - d)Y[d/2 - n 2 + n - m) 
/(nx, n 2 ; n, m) = fRrR . (2.13) 

There are two topologies of two-loop propagator HQET diagrams (figure ||a,6). 
We write the first of them as 
d d ktd d k 2 



-7T 



J (-2w) 2 ^- n3 - n4 -" 5 )J(n 1 , n 2 , n 3 , n 4 , n 5 ) 



n n i n n 2 n n 3 n n4 n n s 
i/j u 2 u 3 u A u 5 

(h+pj-v (k 2 +p)-v 

L>\ — , L) 2 — , 

D 3 = -kf? D 4 = -k 2 2 ? D 5 = -fa - k 2 f . 

(2.14) 

It is symmetric with respect to 1 <-» 3, 2 <-> 4. If one of the indices is zero, it can be 
easily calculated using (|2.12| ) and ( [2.1| ) (figure |5]c,7,e) 

7(m, n 2) n 3 , n 4 , 0) = 7(m, n 3 )I(n 2 , n 4 ) , (2.15) 
7(0, n 2 , n 3 , n 4 , n 5 ) = G(n 3 , n 5 )7(n 2 , n 4 + n 3 + n 5 - d/2) , (2.16) 
J(ni, n 2 , 0, n 4 , n 5 ) = I{ni,n 5 )I(n 2 + n 4 + 2n 5 - n 4 ) (2.17) 

(and symmetric relations). 

Applying the operators d\ ■ {k 4 — k 2 ), d\ ■ k\ and d\ • v to the integrand of ( |2.14j ), 
we obtain [T^ 

[d - ni - n 3 - 2n 5 + + n 3 3+(4- - 5")] 7 = 0, (2.18) 

[d - nx - n 5 - 2n 3 + ml + + n 5 5 + (4~ - 3")] 7 = 0, (2.19) 
[-2nil + + n 3 3+(l" - 1) + n 5 5+(l" - 2 - )] 7 = 0. (2.20) 
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Figure 5: Two-loop HQET propagator diagram 



2 
d 



Of course, more relations are obtained by symmetry. Applying cu-£j and using ho- 
mogeneity in u, we obtain 



\2(d — n 3 — n 4 — n 5 ) — ri\ — n 2 + nxl + + n 2 2 + ] 7 = 



(2.21] 



which is nothing but the sum of ( |2.19|) and its mirror-symmetric. Subtracting the 
2 shifted version of ( p. 21 ) from (|2.18|) , we obtain the most useful relation [16 



d — rii — n 2 — n 3 — 2n 5 + 1 — 
- (2(d - n 3 - n 4 - n 5 ) -n l -n 2 + 1)2" + n 3 3 + (4" - 5") 1 = 0, (2.22) 

which lowers n 2 , n^, and does not raise heavy-quark indices. 

If indices of two adjacent lines are non-positive integers, the integral contains a 
no-scale vacuum subdiagram and hence vanishes. The cases with zero indices are 
given by ( |2.15| ), ( |2.16| ) and ( |2.17| ). When n 2 < 0, it can be raised by ( [2.22| ); the case 
rii < is symmetric. Similarly, if n 3 < 0, it can be raised by ( |2.22| ); the case n 4 < 
is symmetric. When n 5 < and n 3 ^ 1, n 5 can be raised by ( |2. 22| ) (the case n 4 ^ 1 
is symmetric); when n 5 < and n\ 1, 715 can be raised by ( f2.19| ) (the case n 2 ^ 1 
is symmetric). The case < 0, n\ = n\ = n 3 = n 4 = 1 is handled by 



(2.23) 



[(d - 2n b - 4)5+ - 2(d - n 5 - 3)] 1(1, 1, 1, 1, n 5 ) = 

= [(2d - 2n 5 - 7):r 5+ - 3 4+5 + 1~3 + ] 1(1, 1, 1, 1, n 5 ) 



which follows from ( |2.18| ), ( |2.19| ) and ( |2.20| ) at rii = n 2 = n 3 = n 4 = 1 (note that the 
terms on the right-hand side of ( |2.23| ) are trivial for any 715; for n.5 < 0, they vanish). 

We are left with the most important situation when all the indices are positive. 
Applying ( |2.22| ), we reduce n 2 , 114, n 5 until one of them vanishes. Then ( |2.15| ), ( j2.1(j| ) 
and ( |2.17| ) apply. If max(ni, n 3 ) < max(n 2 , n 4 ), it is more efficient to lower rii, n 3 , n 5 . 

In the second topology (figure |5]&), three heavy-quark denominators depend on 
only two variables k\ 2 • v, hence they are linearly dependent. Therefore, there is one 
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scalar product which cannot be expressed via the denominators. Let's define the 
integral 

= - 7rfi (- 2w ) 2(d+n °" n4 " n5)j ^' n2 ' n3 ' n4 ' n5 ' n °) > 
(ki+p)-v (k 2 +p)-v (k 1 + k 2 +p)-v 

~ U! ' ~~ LU ' ~~ U! ' 

L> 4 = -k\ , D s = -k\ , N = 2ki ■ k 2 (2.24) 

(it is symmetric with respect to 1 2, 4 <-» 5). Noting that D\ + D 2 — D 3 = 1, we 



/ 



immediately have [16 



(l-l--2- + 3 _ )J=0. (2.25) 
Applying d\ • ki, d\ • k 2 and d% • v to the integrand, we have 

[d + n - Tii - n 3 - 2n 4 + nil + + n 3 3 + 2~] J = (2.26) 
[m - n 3 - + n 3 3+l" + n 4 4 + 0+ - 2noO~5 - ] J = (2.27) 

[-2^1+ - 2n 3 3 + + n 4 4 + (l~ - 1) + n 0"(2" - 1)] J = . (2.28) 

Homogeneity in u> gives [2(<i+n — n 4 — n 5 )— ni— n3 + n 1 l + + n 2 2 + + ?2 3 3 + ] J = 0, 
which is nothing but the sum of ( |2.26|) and its mirror-symmetric. The boundary 
values of the integral (figure |5|c,e) are 

J(n u n 2 , 0, n 4 , n 5 ; n ) = (5~ - 3 - 4") n °/(ni, n 2 , n 4 , n 5 , 0) (2.29) 
J (0, n 2 , n 3 , n 4 , n 5 ; n ) = {4T - Z + 5~) n °/(n 3 , n 2 , 0, n 5 , n 4 ) (2.30) 

and the symmetric relation for n 2 = 0. If no = 0, 

J (ni,n 2 , 0, n 4 , n 5 ) = /(ni, n 4 )/(n 2 , n 5 ) (2.31) 
J (0, n 2 , n 3 , n 4 , n 5 ) = I(n 3 , n A )I(n 2 + n 3 + 2n 4 - d, n 5 ) . (2.32) 

If n 4 < 0, or 71$ < 0, or two adjacent heavy-quark indices are non-positive, the 
integral vanishes. If any of n\, n 2 , n 3 is negative, it can be raised by fl2.25| ). If all 
of them are positive, we use Q2.25D to lower rii, n 2 or n 3 , until one of these indices 
vanish. 

Instead of using ( |2.29| ), (|2.30| ) when n > 0, we could proceed in another way. If 
n 4 > 1, we lower both n and n 4 using (|2.27|) (the case n 5 > 1 is symmetric). We are 
left with J(rii, n 2 , 0, 1, 1; n ) and J(0, n 2 , n 3 , 1, 1; n ). In the first case, if n\ > 1, we 
lower it using (|2.28|) (the case n 2 > 1 is symmetric). In the second case, if n 3 > 1, 
we lower it using ( [2.28j ); if n 2 > 1, we lower it using the difference of ( |2.28| ) and 
its mirror-symmetric. The two remaining cases (figure |5]c,e) are easily evaluated 
using ( [2.13j ): 

J(l, 1, 0, 1, 1; n ) = (-l) no /(l, 1; n , 0)1(1, 1) , 



n ° ( lYn ! 

J(0, 1, 1, 1, 1; n ) = 1(1, 1; n , 0) £ i--^-J(4 - 2n + / - d, 1 



1=0 
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a b 
Figure 6: Basis two-loop HQET propagator integrals 




All one-loop integrals (figure [|) with integer n^ 2 are proportional to I\ = 1(1, 1), 
the coefficient being a rational function of d. All two-loop integrals with integer 
indices reduce to If (figure |G]a) and I 2 = 1(0,1,1,0,1) (figure ^b), with rational 
coefficients. Here 

In = j. ,\ 9 » r(l + 2ne)T n (l - e) . (2.33) 

(1 - n(d - 2)) 2n 



2.3 Three-loop HQET diagrams with lower-loop insertions 

Diagrams with a two-loop propagator insertion, or with two one-loop insertions, fig- 
ure |7|, are trivially calculated by multiplying the relevant insertionfs] by the one-loop 
HQET integral with non-integer indices, whose values are obvious by dimensionality. 

Next we consider the diagrams obtained from the two-loop HQET diagram of 
figure ||a by adding a single one-loop propagator insertion (figure ||). They are equal 
to the product of the corresponding one-loop integral and the integral of figure |5|a 
with one non-integer index. All relations derived for the diagram of figure [|a are 
valid; however, the non-integer index changes the strategy of their application. 

Let's consider the case of figure [|]a with a non-integer n 3 . When n 5 < 0, it can 
always be raised by (|2.22j ). When n 2 < 0, it can be raised by Q2.22|) , too; when 
rii < 0, it can be raised by the mirror-symmetric relation. When n 4 < 0, it can 
always be raised by ( [2.22|) , too. Finally, when all indices rii, 714, n§ are positive, 
we use ( |2.22| ) to lower n 2 , n 4 or n 5 , until (|2.15| ), Q2.16Q and (|2.17| ) are reached. 

In the case of figure ^b, rii is non-integer. When n 5 < and n 3 7^ 1, we can raise 
n 5 by ( |2.22| ) (the case n 4 7^ 1 is symmetric); when n 5 < and n 2 7^ 1, we can lower 
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a b c 

Figure 9: The diagram of figure |B|5 with a single one-loop insertion 



n 2 by the relation symmetric to (2.19); when < and n 2 = n 3 
use the relation 



n 4 = 1, we can 



[2(d 



n 5 



3)-4+-n 5 5 + l-] J(n 1 ,l,l,l,n 5 ) = 



(2.34) 



which follows from ( |2.20| ) and ( |2.19|) , and then the term with 4 + can be treated as 
above. When n 2 < 0, it can be raised by ( f2.22| ). When n 3 < and n 4 7^ 1, we 
can raise n3 by the relation symmetric to (|2.22|) ; when n 3 < and 7^ 1, we can 
raise n3 by ( |2.19|) ; when n 3 < and n 2 7^ 1, we lower n 2 by ( [2.21 
n 3 < and n 2 = n 4 = n 5 = 1, we use ( |2.22| ) to get a trivial term with n 2 
n 4 < and n 3 7^ 1, we can raise n 4 by (|2.22|); when n 4 < and 7^ 1, we can 



finally, when 
0. When 



raise n 4 by ( |2.19D ; when n 4 < and n 2 7^ 1, we lower n 2 by ( |2.21| ); finally, when 
n 4 < and n 2 = n 3 = = 1, we use ( 2.20| ) to raise n 3 or n^, and proceed as above. 
Finally, when all indices n 2 , n 3 , n 4 , n 5 are positive, we use (|2.22| ) to lower n 2 , n 4 or 
n 5 , until ( |2.15| ), ( |2.16| ) and (|2.17|) are reached. 

In the case of figure |8|c, is non-integer. When n 2 < 0, it can be raised by ( |2.22| ); 
the case n\ < is symmetric. When n 3 < 0, it can be raised by (|2.19|) ; the case 
n 4 < is symmetric. When rix > 1, it can be lowered by (|2.19|) ; the case n 2 > 1 is 



symmetric. When n 3 > 1, it can be lowered by ( 2.20 ), with a trivial additional term 
having n\ = 0; the case n 4 > 1 is symmetric. We are left with n\ = n 2 = n 3 = n 4 = 1; 
the relation (|2.23|) can be used to lower or raise ns, with trivial additional terms. 
An integral with some specific value of 715 (of the form integer plus e) has to be 
considered as a new basis element. This integral has been calculated, exactly in d 
dimensions, in |TU[ in terms of 3F2 hypergeometric functions, using the Hegenbauer 
polynomial technique in the coordinate space |2IJ . 

Finally, we consider the diagrams obtained from the two-loop HQET diagram of 
figure by adding a single one- loop propagator insertion (figure §). The diagram 
of figure Oa, with a non-integer n 4 is calculated exactly as the two-loop one. 
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In the case of figure ^b, n\ is non-integer. When ^3 < or n 2 < 0, they can be 
raised by (|2.25|) . In the case of positive indices, we use 



[d - 7i 3 - 2n A + n 1 l + (2' - 3~) + n 3 3 + 2~] J = 



(2.35) 



(which follows from (2.26) and ( |2.25D ). It either lowers n 2 + n 3 , or, at a fixed n 2 + n 3 , 
lowers n 2 - Therefore, sooner or later, we reach ( 2.29| ) and (|2.30 ). 

In the case of figure |9|c, n 3 is non-integer. If m < or n 2 < 0, they can be 
raised by (|2.25|) . When n 4 > 1, we can lower it or rii by (|2.28|) ; the case n 5 > 1 is 
symmetric. When n\ > 1, it can be lowered by (|2.26|) ; the case n 2 > 1 is symmetric. 
We are left with n\ — n 2 — — 715 — 1; the relation ( |2.25|) can be used to lower or 
raise n 3 , with trivial additional terms. An integral with some specific value of n 3 (of 
the form integer plus 2e) has to be considered as a new basis element. 

It is not difficult to calculate J(l, 1, n 3 , n 4 , n 5 ) for arbitrary n 3 ^ 5 (not necessarily 
integer) in the coordinate space: 

r(n - 2(d -n 1 -n 2 - l))r(d/2 - ni)r(d/2 



J(l, l,n,n h n 2 ) 



n 2 



j ^2(d— ni — n-2)— n— 1 



r(n)r(m)r(n 2 ) 



dt 2 t 2 2 ni ~ 



j. 



\2ri2—d 



(*2 " ^ 



\n-l 



0<ti<i 2 <t 
■ 3-^2 



1, d — 2n 2 , 2ni + n + 1 — d 
n + 1, 2ni+n + 2-rf 



n(2ni + n+l-d)^\ n + 1, 2n x + n + 2 - d J '^' 36 ' ) 

All diagrams considered in this section are particular cases of the generic three- 
loop topologies, which will be discussed in section ^|, when some lines are shrunk 
(i.e., some indices vanish). Therefore, we don't consider diagrams with numerators 
here: numerators should be dealt with in the context of generic topologies, and 
the formulae of this section are used only as boundary values for the corresponding 
recurrence relations, after elimination of numerators. 

All three-loop HQET propagator integrals with lower-loop propagator insertions 
are linear combinations of 7 basis integrals (figure |T0|a-ff), coefficients being rational 
functions of d. The basis integrals of figure [TTja (If), figure |T0| fe (hh), figure |T0| c (7 3 ), 
figure [[(](/ (I 3 If/I 2 ), and figure [TTJe (I 3 G\/G 2 ) are known exactly in d dimensions in 
terms of V functions. Those of figure \£(\f,g contain hypergeometric 3 F 2 functions of 
the unit argument. Several terms of their expansion in e can be obtained using the 
methods which were recently developed in . As we shall see in section |3|, there is 
only one additional basis integral, figure |T01/l 



3. Proper three-loop HQET propagator diagrams 

In the massless case 0, there are only 3 topologies of proper three-loop propagator 
diagrams: Mercedez, Ladder, and Non-planar (plus their reduced forms obtained 
by shrinking some lines). Now we have 10 topologies instead (figure [TT| ). Each of 
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f 9 h 

Figure 10: Basis tree-loop HQET propagator integrals 

them has 8 propagators. In the diagrams with 4 heavy-quark lines (figure |TT|e- 
g), there is one linear dependence between their denominators; with 5 heavy-quark 
lines (figure [ll]/i-j) - - 2 dependences. There are 9 independent scalar products 
of 3 loop momenta &i,2,3 and the 4-velocity v. Therefore, in the diagrams with 
two or three heavy-quark lines (figure [TT|a— d) , there is one scalar product in the 
numerator which cannot be cancelled against the denominators; with 4 heavy-quark 
lines (figure [TTje — ) — two scalar products; with 5 heavy-quark lines (figure P]/t— j) 
— three scalar products. 

When calculating these diagrams using recurrence relations, some indices may 
vanish. This corresponds to shrinking the corresponding lines. In some cases, this 
results in diagrams with lower-loop propagator insertions, which were calculated in 
section [2]. The diagrams of figure [T2] are still non-trivial. 

3.1 Diagram with two heavy-quark lines 

Let's consider the diagram of figure [TT]a first. We define 
N no d d k 1 d d k 2 d d k 3 



n n i n™2 n"s n n 4 n™5 n n e n™ 7 n n » 

U X U 2 LJ^ J_7 4 J_7 5 ±7g Urj U % 

" I a {ni, n 2 , n 3 , n 4 , n 5 , n 6 , n 7 , n 8 ; n ) 









N = 


k 3 ■ v 


Di 






00 


D 3 = 


-k 2 




D 6 = 


ft 3 J 


D 7 



(h+p)-v (k 2 +p)-v 

; U 2 — , 

LO 00 

-k 2 2 , D 5 = —(hi - k 2 ) 2 , 

-(k 3 + h) 2 , D 8 = -(k 3 + k 2 ) 2 . (3.1) 

This integral is mirror-symmetric with respect to 1 «-> 2, 3 <-> 4, 7 <-> 8. It vanishes 
when the indices of the following groups of lines are non-positive: 12, 67, 68, 78, 375, 
485, 315, 425, 364, 137, 248, 157, or 258. 
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12 1 3 2 

a b 




13 2 13 2 

c d 




h % 3 

Figure 11: Topologies of proper three- loop HQET propagator diagrams: Mercedez (a, b, 
f, g, i), Ladder (c, e, h), Non-planar (d, j) 




d e f 

Figure 12: Topologies of non-trivial three-loop diagrams with a shrunk line 



12 



First we are going to get rid of the numerator. When no > and n 7 ^ 1, we can 
lower no by 

[-272x1+ + n 3 S + (l- - 1) + n 5 5 + (l- - 2") + n 7 7+(l~ - 1 + 0+)] 4 = 0, (3.2) 



which is obtained by applying di -v to the integrand of ( |3.1| ); the case n > 0, n 8 ^ 1 
is symmetric. When n > and n 6 7^ 1, we can lower no by 



2 [3d + n - ni - n 2 - 2(n 3 + n 4 + n 5 + n 6 + n 7 + n 8 )] + 



la = , (3.3) 



+n 3 3 + (l~ - 1) + n 4 4+(2~ - 1) - n 6 6 + + - 2n 0~ 
which is the {d\ + d 2 — 83) ■ v relation simplified using the homogeneity relation 

[3d + n - rii - n 2 - 2(n 3 + n 4 + n 5 + n 6 + n 7 + n 8 ) + nil + + n 2 2 + ] J a = . (3.4) 
When n > and n\ ^ 1, we can lower n by 

d — ni — n 3 — n 5 — 2n 7 + nil + (l — + ) + 



+n 3 3 + (6~ - 7-) + n 5 5 + (8" 







(3.5) 



which is obtained by applying d\ ■ (k\ + k 3 ) to the integrand of ( |3.1| ); the case n > 0, 
n 2 7^ 1 is symmetric. We are left with n > 0, n 4 = n 2 = n 6 = n 7 = n 8 = 1; we use 
$3 ■ (^3 + k%) relation 



d + n - n 6 - n 8 - 2n 7 + n (1 - 1) + 



+n 6 6 + (3--7-)+n 8 8 + (5" 







(3.6) 



to lower rio or raise n 6 or n 8 , and apply the method described above again. 

Now we shall discuss the integral I a {ni, n 2 , n 3 , n 4 , n 5 , n 6 , n 7 , n 8 ) without numer- 
ator (n = 0). Applying <9 3 • fc 3 , <9i • k±, d\ ■ {k\ — k 2 ), (<9i + d 2 — <9 3 ) ■ k\ to the integrand 
of ( |3.1| ), we obtain the recurrence relations 

[d - n 7 - n 8 - 2n 6 + n 7 7 + (3~ - 6") + n 8 8 + (4~ - 6")] 4 = 0, 
d — rii — n 5 — n 7 — 2n 3 + nil + + 



+ n 5 5 + (4~ - 3") + n 7 7 + (6~ - 3" 
d — rii — n 3 — n 7 — 2n 5 + n^+2 - - 
+ n 3 3+(4~ - 5") + n 7 7+(8~ - 5 
2(d — ns — n 7 — n 8 ) — n 2 — n 4 — n 6 + n 2 2 + l~ - 
+ n 4 4 + (3" - 5") + n 6 6 + (3~ - T 
where the last relation was simplified using 







(3.7) 
(3.8) 



4 = 



(3.9) 



(3.10) 



13 



The cases n% = 0, n 2 = 0, tiq = 0, n 7 = 0, n§ = are trivial. When n x < 0, it 
can be raised by 

[3d - ni — n 2 — 2(n 3 + n 4 + n 5 + n 6 + n 7 + n 8 )] / a = (3-H) 
= [d - m - n 2 - n 4 - n 8 - 2n 5 + n 4 4+(3~ - 5") + n 8 8 + (7~ - 5")] 1+4 , 



which is the difference of (3.4) and 1 + shifted version of the relation symmetric 
to (|3.9| ); the case n 2 < is symmetric. When n 6 < and n 7 7^ 1, we can raise 
Uq by (|3.8| ) (the case n 8 7^ 1 is symmetric); when n 6 < and n 7 = n 8 = 1, we 
can raise n.6 or n% by ( |3.6| ). When 727 < and n% 7^ 1, we can raise n 7 by the 
relation symmetric to ( |3.9|) ; when n 7 < and n 6 7^ 1, we can raise n 7 by ( |3.10| ); 
when n 7 < and n 6 = n 8 = 1, we can raise n 7 or n 6 by the relation sym- 
metric to (|3.6| ). The case n§ < is symmetric. The case n 3 = (figure 12a. 



J a (ni, n 2 , n 4 , n^, n 6 , n 7 , n 8 )) will be considered later in this section; the case n 4 = 
is symmetric. When n 3 < and n 6 7^ 1, we can raise n 3 by (|3.6|) ; when n 3 < 
and n 7 7^ 1, we can raise n 3 by ( |3.7|) ; when n 3 < and n 6 = n 7 = 1, we can 
raise or n 7 by the relation symmetric to ( |3.6|) . The case n 4 < is symmet- 
ric. The case n 5 = (figure |l"2"| fe, Jb(n x , n 2 , n 3 , n 4 , n 6 , n 7 , n 8 )) will be considered 
later in this section. When n 5 < and n 8 7^ 1, we can raise n 5 by ( |3.6| ) (the 
case n 7 7^ 1 is symmetric); when < and n 7 = rig = 1, we can raise n 7 or 
n 8 by (|3.7|). When all the indices are positive, we can use ( |3.7| ) to kill one of the 
lines 3, 4, 6; or we can use Q3.6p to kill one of the lines 3, 5, 7; or we can use 
the symmetric relation to kill one of the lines 4, 5, 8; or we can use ( |3.10|) to kill 
one of the lines 1, 3, 5, 7; or we can use the symmetric relation to kill one of the 
lines 2, 4, 5, 8. 

Now we consider J a {n x , n 2 , n 4 , n 5 , n 6 , n 7 , n 8 ) = I a {n x , n 2 , 0, n 3 , n 4 , n 5 , n 6 , n 7 , n 8 ) 
(figure [12]a). This integral vanishes when the indices of the following groups of lines 
are non-positive: 12, 67, 68, 78, 57, 15, 17, 46, 485, 425, 248, or 258. If any index 
is zero, the integral becomes trivial. If n x < 0, n 2 < 0, n e < 0, n 7 < 0, or n s < 0, 
we just consider J a as I a with n 3 = and proceed as usual; the integral reduces to 
trivial ones not including J a . When n 4 > 0, we can kill one of the lines 2, 4, 8 using 
the relation symmetric to (|3.10|) . When n 5 > 0, we can kill one of the lines 2, 5, 8 
using (|3~9| ). When both n 4 and are negative, we proceed as follows. If n 6 7^ 1, 
n 4 can be raised by the relation symmetric to ( |3.10|) ; if n 7 7^ 1, n 5 can be raised 
by (|3.9|) ; n 2 can be lowered down to 1 using (|3.4|) . The remaining integral can be 
calculated by a repeated use of ( j2.13j ): 

J a (ni, 1, -|n 4 |, -\n 5 \, 1, l,n 8 ) = 





n 4 


l 


n 5 


i 


(l«5 


-hy.mlih -2m)!(|ra 4 
,m)I(l,n 8 - |n 4 - |n 5 


| -Z 2 )!Z 3 !(/ 2 -Z3)! 
+ h + 1 2 — m; l 2l 



E 

h,l 2 ,l3,m 



x J(-2(d + |n 4 | + |n 5 | - ra 8 ) + n x + Z 3 + 3, 1) 



x 
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Finally, we consider J&(ni, n 2 , n 3 , 714, n 6 , n 7 , n 8 ) = I a (n 1 ,n 2 ,n 3 , 71^,0, n 6 ,n 7 ,n 8 ) 
(figure |l2|fe). This integral is mirror-symmetric; it vanishes when the indices of the 
following groups of lines are non-positive: 12, 67, 68, 78, 37, 48, 13, 17, 24, 28, or 
364. If any index is zero, the integral becomes trivial. If n\ < 0, n 2 < 0, n 3 < 0, 
n 4 < 0, tiq < 0, n 7 < 0, or n 8 < 0, we just consider as I a with n 5 = and proceed 
as usual; the integral reduces to trivial ones not including J b (though, possibly, 
including J a ). When all the indices are positive, one of the lines 3, 4, 6 can be killed 

by Q. 



3.2 Mercedez with three heavy-quark lines 



Now let's consider the diagram of figure [Tip. We define 
N^dtkxdPkifPkz 



j D™ 1 DV? D% 3 1?4 4 5 Dq 6 D 7 7 D n& 



r 3d/2. 



8 

8 



-2uj) 3d+2no - 2 ^=i n ' I, 



(k!+p)-v (k 2 +p)-v _(k 3 +p)-v 

U \ — , L> 2 — , U 3 — , 

Di = -k\ r D 5 = -k\ ^ D 6 = -(A* - hf , 

D 7 = -(k 2 - hf , D 8 = -fa - k 2 f , N = -kl . 

(3.12) 

This integral is mirror-symmetric with respect to 1 <-> 2, 4 <->• 5, 6 <->• 7. It vanishes 
when the indices of the following groups of lines are non-positive: 36, 37, 67, 123, 
148, 168, 416, 486, 258, 278, 527, or 587. 

First we are going to get rid of the numerator. When n& ^ 1, we can lower n 
by the di ■ (k% — k 3 ) relation 

[d - nx - n 4 - n s - 2n 6 + ml + 3~ + n 4 4+(0 + - 6") + n 8 8 + (7~ - 6")] I b = ; 

(3.13) 

the case n 5 ^ is symmetric. When n% ^ 1, we can lower n by the di ■ k% relation 

[d - rix - n 6 - n 8 - 2n 4 + n x l + + n 6 6 + (0 + - 4~) + n 8 8+(5~ - 4")] I b = 0; 

(3.14) 

the case n 7 ^ 1 is symmetric. When 713 7^ 1, we can lower no or raise or 77,7 by the 
d 3 ■ v relation 

[-2n 3 3 + + noO~(l - 3") + n 6 6 + (3~ - 1") + n 7 7+(3~ - 2")] I b = . (3.15) 

Finally, we can lower or raise n 3 or n 7 by the d% ■ (k 3 — k\) relation 

[d + n - n 3 - n 7 - 2n 6 + n 3 3 + l" + n (T(6~ - 4 _ ) + n 7 7 + (8" - 6")] I b = . 

(3.16) 
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Now we shall discuss the integral 7&(ni, n 2 , n 3 , n 4 , n.5, n$, n?, n 8 ) without numer- 
ator (no = 0). Applying d\ • v, d± ■ {k\ — k 2 ), (d\ + d 2 + <9 3 ) ■ k\ to the integrand 
of ( |3.12| ), we obtain the recurrence relations 

[-2n a l + + n 4 4 + (l- - 1) + n 6 6 + (l- - 3~) + n 8 8 + (l~ - 2~)] I b = , (3.17) 

d — ni — — riQ — 2ns + n x l + 2~ + 



+ n 4 4+(5~ - 8") + n 6 6 + (7- - 8") 
2(d — n e — n 7 — n 8 ) - n 2 - n 3 - n 5 + 



+ (n 2 2 + +n 3 3 + )l- +n 5 5 + (4" 



8" 



0, (3.18) 



0, (3.19) 



where the last relation was simplified using the homogeneity relation 



3d — n 1 — n 2 — n 3 — 2(n 4 + n 5 + n 6 + n 7 + n 8 ) + 

+ nil + + n 2 2 + + n 3 3 + 



h = 0. 



(3.20) 



Using it to simplify the sum of ( |3.17| ), its symmetric and ( |3.15| ), we get 



2 [3d — n x — n 2 — n 3 — 2(n 4 + n 5 + n 6 + n 7 + n$)] + 
+ n 4 4+(l- - 1) + n 5 5 + (2- - 1) 



h = 



(3.21) 



The cases n 3 = 0, n 4 = 0, ns = 0, n 6 = 0, 717 = are trivial. When n 3 < 
and «6 7^ 1, we can raise n 3 by (|3.17D (the case n 7 7^ 1 is symmetric); when n 3 < 
and tie — n 7 — 1, we can raise n 3 or 77,7 by (|3~T6| ). When n A < and n 5 7^ 1, 
we can raise n 4 by the relation symmetric to (|3.18| ); when n 4 < and ns — 1, we 
can raise n 4 or n 5 by ( |3.21|) . The case n 5 < is symmetric. When n 6 < and 
n 7 7^ 1, we can raise n§ by ( |3.16| ); when n 6 < and n 3 7^ 1, n 7 = 1, we can raise 
n6 or n 7 by ( |3.15|) ; when n@ < and n 3 = n7 = 1, we can raise n§ or n 3 by the 
relation symmetric to ( |3.16| ). The case n7 < is symmetric. The case n\ = 

is symmetric. When 



is J a (n 3 ,n 2} n 5 ,n 7} n 4} n 6 ,n s ) (figure [TJa, section pTT|), n 2 

< and n 8 7^ 1, we can raise n x by the relation symmetric to ( |3.17| ); when nx < 
and n6 7^ 1, we can raise n\ by fl3.15|) ; when n\ < and n 4 7^ 1, we can raise rii 
by ( |3.21| ); when n\ < and n 2 7^ 1, we can raise n\ by the relation symmetric 
to ( |3.18|) ; when n x < and n 3 7^ 1, we can raise n\ by ( |3.16|) ; when ni < and 
n 2 = n 3 = n 4 = ne = n$ = 1, we can raise ni, n 2 or n 3 by ( |3.20 ). The case 
n 2 < is symmetric. The case n 8 = (figure [I2|c, J c (ni,n 2 , n 3 , n 4 , n 5 , n 6 , n 7 )) will be 



considered later in this section. When n 8 < and % 7^ 1, we can raise n 8 by ( |3.16| ) 
(the case n 6 7^ 1 is symmetric); when n s < and n 5 7^ 1, we can raise n 8 by ([ 
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(the case n 4 7^ 1 is symmetric); when rig < and n\ 7^ 1, we can raise ng, n 4 or 
n6 by (|3.17|) (the case n 2 7^ 1 is symmetric); when ng < and n 3 7^ 1, we can raise 
n 6 or n 7 by ( |3.15| ); when n 8 < and n x = n 2 = n 3 = n 4 = n 5 = n 6 = n 7 = 1, 
we can raise ni, n2 or n 3 by ( |3.20| ). When all the indices are positive, we can kill 
one of the lines 1, 6, 8 by ( |3.16| ), or one of the lines 2, 7, 8 by its mirror-symmetric 
relation. 

Now we consider J c (n 1 ,n 2 , n 3 , n 4 , n 5 , n 6 , n 7 ) = I b (n l} n 2} n 3 , n 4 , n 5 , n 6 , n 7 , n 8 ) (fig- 
ure 0c). It is mirror-symmetric; it vanishes if any two indices of n 4 , n 5 , n 6 , n 7 are 
non-positive, or n 1; n 2 , n 3 are all non-positive. If any index is zero, the integral 
becomes trivial. If any index is negative, we just consider J c as I b with n$ = and 
proceed as usual. Using the d\ ■ v and <9 3 • v relations 



[-2^1+ + n 4 4 + (l- - 1) + n 6 6+(l" - 3")] J c = 
[-2n 3 3 + + n 6 6 + (3" - 1") + n 7 7 + (3" - 2")] J c = 



(3.22) 
(3.23) 



we can lower ni, n 2 (symmetric case) and n 3 down to 1. 

We are left with J c (l, 1, 1, n 4 , n 5 , n 6 , n 7 ). It is rather difficult to apply the stan- 
dard techniques to this integral, because each operator c\ • kj produces a scalar 
product which cannot be expressed via the denominators. For example, d\ ■ k\ gives 
the term — j^2ki ■ k 3 , and d\ ■ (ki — k 3 ) gives the term — jf2k\ ■ k 3 . We can cancel 
2ki ■ k 3 by forming the difference 



(n 4 - 1)6-0! • h - (n 6 - l)4-d 1 ■ (h - k 3 ) 



which results in 



n 4 — 1) (d — rii — 2n 4 + nxl + ) 6 — (n 6 — l)(d — n\ — 2n 6 + riil + 3 ) 4 

Another useful combination is 

(n 4 - 1)5- (d 1 + d 2 + d 3 ) ■ k x - (n 5 - 1)4- (d 1 + d 2 + d 3 ) ■ k 2 



J c = 0. 
(3.24) 



(it can be simplified by the homogeneity relation). Using (|3.22j ) and (|3.23| ), we obtain 
at rii = n 2 = n 3 = 1 

(n 4 - 1) [2(d - 2n 4 - 1) - n 4 4 + (l - 1~)] 6" - 
2(n 6 - - 2n 6 - 1 + l + 3-)4- + (n 4 - l)(n 6 - 1)(1" - 3 _ )1 J c = , (3.25) 

(n 4 - 1) [2(d - n 5 - n 6 - n 7 ) + (2+ + 3 + )l-] 5~ - 
- (n 5 - 1) [2(d - n 4 - n 6 - n 7 ) + (1+ + 3 + )2"] 4~1 J c = . (3.26) 
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One more useful relation is obtained by adding ([3.22 ), its mirror-symmetric, for- 
mula (|3.23|) and using the homogeneity relation: 



2 [3d — ni — n 2 — n 3 — 2(n 4 + n 5 + n 6 + n.7)] + 
+ n 4 4+(l- - 1) +n 5 5 + (2- - 1) 



Jr = 0. 



(3.27) 



When Tie > 1, we can lower it down to 1 (raising n 4 ) by (|3.25 ); the case n 7 > 1 is 
symmetric. Finally, we can lower n 4 and 715 down to 1 by ( |3.26| ) together with ( 3.27Q . 

The integral J c (l, 1, 1, 1, 1, 1, 1) cannot be reduced to simpler ones (in contrast 
to the massless case [|J), and should be considered as a new basis integral. Its value 
is currently unknown, and its calculation is highly non-trivial. 

3.3 Ladder with three heavy-quark lines 

Now let's consider the diagram of figure |TT] c. We define 

N no d d k 1 d d k 2 d d k 3 
nni n na n™ 3 n" 4 n™s r> n s n" 7 n n & ~~ 

u l u 2 u 3 u A u 5 lj 6 u 7 u % 

= -nr 3d/2 (-2uj) 3d+2no - 2 ^ n *I c (n u n 2 , n 3 , ra 4 , n 5 , n 6 , n 7 , n 8 ; no) , 



D 7 



(/ci +p) • u 



a; 



-(h~k 3 ) 2 



D 2 



(k 2 + p) ■ v 



-k 2 
-k 2 

K 3 ) 



^3 



(k 3 +p) ■ v 



to 

-(h - h- 

N = 2ki ■ k 2 



2 

3; , 



(3.28) 

This integral is mirror-symmetric with respect to 1 <-» 2, 4 <-» 5, 6 <-» 7. It vanishes 
when the indices of the following groups of lines are non-positive: 14, 16, 46, 25, 27, 
57, 132, 637, 368, 378, 687, 485, 487, 586. 

First we are going to get rid of the numerator. When n 7 7^ 1, we can lower no 
by the d 3 ■ (k 3 — ki) relation 



d — n 3 — n-j — n 8 — 2n 6 + n 3 3 + l 4 
+ n 7 7 + (0 + + 4- + 5- - 6-) + n 8 8 + (4- - 6 
the case n 6 7^ is symmetric. When n 5 7^ 1, we can lower n by 
2(d — n 5 — tiq — n 7 ) + n — n 2 — n 3 — n 8 + 2n _ 4~ -f 
+ (n 2 2 + + n 3 3 + )l- - n 5 5 + + + n 8 8 + (4~ - 6" 



0: 



(3.29) 



0. 



(3.30) 



which is the (di + d 2 + d 3 ) • ki relation simplified by the homogeneity relation, the 
case n 4 7^ 1 is symmetric. When n\ 7^ 1, we can lower uq or raise n 4 or by the 
d\ ■ v relation 



[-272x1+ + n Q-(2- - 1) + n 4 4 + (l - - 1) + n 6 6 + (l~ - 3")] I c = 



(3.31) 
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the case n 2 7^ 1 is symmetric. Finally, we can raise n\ or n§ by the d\ ■ k\ relation 



[d + n - ni - n 6 - 2n 4 + 71x1+ + n 6 6 + (8~ - 4")] J c = . (3.32) 

Now we shall discuss the integral J c (ni, n 2 , n 3 , n 4 , n 5 , Uq, n 7 , n 8 ) without numer- 
ator (no = 0). Applying <9 3 ■ v , d\ ■ (k\ — fc 3 ) to the integrand of ( |3.28|) , we obtain the 
recurrence relations 

[-2n 3 3 + + n 6 6 + (3~ - 1") + n 7 7 + (Z- - 2~) + n 8 8+(3" - 1)] I c = , (3.33) 
[d - m - n 4 - 2n 6 + + n 4 4+(8~ - 6")] J c = . (3.34) 

Homogeneity in u gives the relation identical with (|3.20| ). 

The cases n\ = 0, n 2 = 0, n 4 = 0, = 0, n 6 = 0, n 7 = are trivial. When 



ni < and n e 7^ 1, we can raise n\ by ( |3.33| ); when n\ < and n 4 7^ 1, we can raise 
n\ using the sum of ( |3.31| ) and ( |3.33| ); when n\ < and n 4 = ne = 1, we can raise 



n\ or n6 by ( |3.32| ). The case n 2 < is symmetric. When n 4 < and n$ 7^ 1, we can 
raise n 4 by (|3.32|) ; when n 4 < and n\ 7^ 1, — 1, we can raise n 4 or n 6 by (|3.31 



when n 4 < and rii = = 1, we can raise n 4 or ni by (|3.34|) . The case < is 
symmetric. When n@ < and n 4 7^ 1, we can raise n§ by ( |3.34|) ; when uq < and 



ni 7^ 1, n 4 = 1, we can raise n 6 or n 4 by ( |3.31| ); when n 6 < and n\ — n 4 = 1, 
we can raise n 6 or ni by ( |3.32| ). The case n 7 < is symmetric. The case n 3 = 
is Jbijii, n 2 , n 4 , ns, n 8 , ng, n 7 ) (figure |I~2|a, section |3~Tl"| ). When n 3 < and n6 7^ 1, we 
can raise n 3 by ( |3.31| ) (the case n 7 7^ 1 is symmetric); when n 3 < and n 8 7^ 1, we 
can raise n 3 by the relation 

2 [3d - n\ — n 2 — n 3 — 2(n 4 + n 5 + n 6 + n 7 + n 8 )] + 

+ n 4 4 + (l _ - 1) + n 5 5+(2~ - 1) + n 8 8 + (3~ - 1) I c = , (3.35) 

which is ( p.31| ) plus its mirror-symmetric plus (|3.33|) simplified by the homogeneity 
relation; when n 3 < and n\ 7^ 1, we can raise n 3 by ( |3.34| ) (the case n 2 7^ 1 
is symmetric); when n 3 < and ni = n 2 = n 6 = n 8 = 1, we can raise n 1; n 2 
or n 3 by the homogeneity relation. The case n 8 = is J c {ni,n 2 , n 3 ,n 4 , ns,n 6 , n 7 ) 
(figure |T2"|c, section ^.2[ ). When n 8 < and n§ 7^ 1, we can raise n 8 by ( |3.32| ) (the 
case n 7 7^ 1 is symmetric); when n 8 < and n 4 7^ 1, we can raise n 8 by ( p.34|) (the 
case n 5 7^ 1 is symmetric); when n 8 < and n 3 7^ 1, n 4 = n 6 = n 7 = 1, we can 
raise n 8 , n 6 or n 7 by (|3.33|) ; when n 8 < and n\ 7^ 1, n 3 = n 4 = n 5 = n 6 = n 7 = 1, 
we can raise n 4 or by (|3.31|) (the case n 2 7^ 1 is symmetric); when n 8 < and 
m — n 2 — n 3 = n 4 = ns = ne — n 7 = 1, we can raise m, r?-2 or n 3 the homogeneity 
relation. When all the indices are positive, we can kill one of the lines 3, 6, 8 by ( |3.34j ), 
or one of the lines 3, 7, 8 by its mirror-symmetric relation. 
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3.4 Non-planar diagram with three heavy-quark lines 

Now let's consider the diagram of figure |TT]<i. We define 

N no d d kid d k 2 d d k 3 



I 



pj-ni p,U2 pins p,U4, p,ri3 T-)ne ryn nns 

u x u 2 u 3 ±7 4 J_y 5 u 6 u 7 u % 

= -in 3d/2 (-2uj) 3d+2no - 2 ^ n *I d (n u n 2 , n 3 , n 4 , n 5 , n 6 , n 7 , n 8 ; n ) 



Di 



(fa +p) ■ v 



UJ 



D 4 = -k\ , 

D 7 = -(k 2 -k 3 ) 2 



D 2 



(k 2 + p) ■ v 



UJ 



-k 2 

-P 



D 3 



(ki + k 2 - k 3 +p) -v 



-(fa - k 3 
N = 2fa ■ k 2 . 



UJ 
2 



(3.36) 

This integral is mirror-symmetric with respect to 1 <-> 2, 4 <-* 5, 6 <-> 7. It vanishes 
when the indices of the following groups of lines are non-positive: 46, 57, 132, 485, 
487, 586, 687, 314, 136, 325, 237, 148, 178, 258, 268, 417, 526, 637, 368, 378. 
When n\ < 0, n 2 < or n 3 < 0, we use 

Id(-\ni\, n 2 , n 3 , n 4 , n 5 , n 6 , n 7 , n 8 ; n ) = 

= (0+ + l-)l ni l(5~ + 6- - 3 - 8-) no I a (n 3 , n 2 , 0, n 5 , n 6 , n 7 , n 4 , n 8 ) , (3.37) 
Id{ni,n 2 , -\n 3 \,n4,n 5 ,n 6 ,n 7 ,n 8 ;n ) = 

= (0+ + I" + 2 - 1)^(5- - 3 - 4T) n H a (n u n 2 , n 4 , n B , 0, n 8 , n 6 , n 7 ) , (3.38) 



and relation symmetric to ( |3.37| ) . Applying d 2 -v and —d 3 -v to the integrand of ( |3.36| ), 
we obtain 



2n 2 2 + + 2n 3 3 + + n 5 5 + (l - 2~) + n 7 7 + (l _ - 3" 







(3.39) 



2n 3 3 + + n 6 6 + (2- - 3") + n 7 7 + (l" - 3") + 

+ n 8 8 + (l- + 2--3--l)]/ d = 0. (3.40) 

When n\ > 1, we can lower it by 

2[3d -ni-n 2 -n 3 - 2(n 4 + n 5 + n G + n 7 + n 8 - n )] + 

+ n 0-(l- - 1) + 2n 1 l + + n 5 5 + (2- - 1) +n 7 7+(3" - J d = 0, (3.41) 

which is twice the homogeneity relation 

3d - nx - n 2 - n 3 - 2(n 4 + n 5 + n 6 + n 7 + n 8 - n ) + 

+ ni l + + n 2 2 + + n 3 3 + ]l d = (3.42) 



minus ( |3.39| ); the case n 2 > 1 is symmetric. When n 3 > 1, we can lower it by ( |3.40| ). 
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We are left with 1^(1, 1, 1, n 4 , n 5 , Uq, n 7 , n 8 ; n ), and now are going to get rid of 
the numerator. Applying d 3 ■ (k 3 — k 2 ) and d\ • (k\ — k 3 ) to the integrand of (|3.36| ), 
we obtain 



d — n 3 — tiq — n$ — 2n 7 + n 3 3 + l -\ 
+ n 6 6 + (4~ + 5" - 7~ + 0+) + n 8 8 + (5~ - 7~ 

d + n — n 3 — n 4 — 2n 6 + n - (5~ — 7~ + 8~) 4 

+ ni l + (2- - 3") + n 3 3+2~ + n 4 4 + (8~ - 6" 

When n 4 7^ 1, we can lower uq by 

2(d — n 4 — n 6 — n 7 ) + n — n\ — n 2 — n 3 — n 8 + 1 — 
— [3d - ni — n 2 — n 3 + 1 - 2(n 4 + n 5 + n e + n 7 + n$ - n )]2~ + 

+ 2n CT5- - n 4 4+0 + + n 8 8 + (5~ - 7") 



0. 



(3.43) 



(3.44) 



0, (3.45) 



which is the — (<9i + d 2 + d 3 ) ■ k 2 relation plus (|3.42|) minus 2~ shifted ( p.42|) ; the 
case 77,5 7^ 1 is symmetric. When n§ 7^ 1, we can lower n by ( |3.43| ); the case n 7 7^ 1 
is symmetric. Finally, when n 4 = n§ = = n 7 = 1, we can lower no or raise n 4 



by O). 

Now we shall discuss the integral Id(n>i, n 2 , n 3 , n 4 , n 5 , n 6 , n 7 , n 8 ) without numer- 
ator (n = 0). Applying di ■ ki, d 3 ■ k 3 , (<9 2 + d 3 ) ■ k 3 to the integrand of ( [3.361 ), we 
obtain 

d - ni - n 6 - 2n 4 + ml + + n 3 3 + (l - 1~) + n 6 6 + (8~ -4") J d = , (3.46) 

d - n 3 - n 6 - n 7 - 2n s + n 3 3 + (l~ + 2 - 1) + 

+ n 6 6 + (4- - 8-) + n 7 7 + (5- - 8~)1 7 d = , (3.47) 

d - n 2 - n 5 - n 6 - 2n 8 + n 2 2 + (l + 3 - 1 _ ) + 

+ n 5 5+(7- - 8") + n 6 6 + (4- - 8")] J d = . (3.48) 

Some other useful relations are 

2n 2 2 + + n 5 5+(l - 2~) + n 6 6 + (3~ - 2~) + n 8 8+(l + 3 -1- 2")1 J d = , (3.49) 



2(d — n 5 — n 7 — n 8 )—n 2 — n 3 — n 6 + n 2 2 + + n 3 3 + l +n 6 6 + (4 — 8" 



0, (3.50) 



the differences of flOJ) and flOOD , and of (ggg) and ( ggp . 

The case n\ = is J a (n 3 , n 2 , n^, n 6 , n 7 , n 4 , n 8 ) (figure [L2]a, section pTT[ n 2 = is 
symmetric); the case n 3 = is Jb(ni, n 2 , n 4 , ns, n 8 , n 6 , n 7 ) (figure |12|fe, section |3~T|) ; 
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the case n 8 = is J c (ni, n 2 , n 3 , n 4 , n 5 , 77,7, n 6 ) (figure [T^c, section |3?2j). The cases 
n 6 = (Jd(rii,n 3 ,n 2 ,n 4 ,n 7 ,n 5 ,n 8 ), figure |X2|cZ ; 7^7 = is symmetric) and n 4 = 
( J e (ni, n 3 , n 2 , n 6 , n 8 , n 5 , n 7 ), figure [L^e; n 5 = is symmetric) will be discussed later 
in this section. When n 2 < 0, it can be raised using 



[3d - n x - n 2 - n 3 - 2(n 4 + n 5 + n 6 + n 7 + ng)]-^ = 
= [d — n 2 — n 3 — n 4 — 2n G — riil + 3~ + n 4 4 + (8~ — 6 



}2+L 



(3.51) 



which is (|3.42| ) minus 2 + shifted ( [3.44 ); the case n\ < is symmetric. When n 3 < 
and Hi 7^ 1, we can raise n 3 by ( |3.44|) (n 2 7^ 1 is symmetric); when n 3 < and 
n 7 7^ 1, we can raise n 3 by ( |3.39| ) (?^6 7^ 1 is symmetric); when n 3 < and n 8 7^ 1, 
we can raise n 3 using 



2[3d — ni — n 2 — n 3 — 2(n A + n 5 + n 6 + n 7 + n 8 )]+ 



+ n 4 4+(l- - 1) +n 5 5+(2- 



n 8 8 + (r 



0, (3.52) 



which is ( p.39| ) plus its symmetric minus ( |3.4(J| ) simplified by ( |3.42| ); when n 3 < 
and ni = n 2 = n§ = n 7 = n 8 = 1, we can raise n 3 , ri\ or n 2 by ( |3.42| ). When ri\ > 1, 
it can be lowered by (|3.41| ); the case n 2 > 1 is symmetric. When n 3 > 1, it can be 
lowered by (|3.40| ). 

We are left with 7^(1, 1, 1, n 4 , n 5 , n 6 , n 7 , n 8 ). When n 4 < 0, it can be raised 
by (|3.44j ); the case < is symmetric. When tiq < and n 4 7^ 1, we can raise uq 
by ( |3.44|) ; when n% < and n 4 = 1, we can raise n§ or n 4 using 

2{d - n x - n 6 - 2n 4 ) - 2n 3 3 + l" + 7i 4 4 + (l~ - 1) + 



+ n 6 6 + [2(8~ -4") +3" -2 



(3.53) 



which is twice (|3.46| ) minus the relation symmetric to ( |3.39| ). The case n 7 < is 
symmetric. When n 8 < and n 4 7^ 1, we can raise n 8 by Q3.44|) (77,5 7^ 1 is symmetric); 
when n 8 < and ri6 7^ 1 ; we can raise n 8 by ( p. 53 ) (n 7 7^ 1 is symmetric); when 
n 8 < and n\ = n 2 = n 3 = n 4 = n 5 = n G = n 7 = 1, we can raise n 8 , n 4 or n 5 
by ( |3.52|) . When all the indices are positive, we can kill the line 6 or 8 using (|3.44| ). 

Thus, the non-planar diagram reduces to planar ones, in contrast to the massless 
case where such a reduction is impossible || |22| . 

Now we consider J d {n u n 2 , n 3 , n 4 , n 5 , n 7 , n 8 ) = I d (rii, n 2 , n 3 , n 4 , n 5 , 0, n 7 , n 8 ) (fig- 
ure ^d). This integral vanishes when indices of the following groups of lines are 
non-positive: 4, 13, 23, 57, 58, 78, 37, 38, 25, 28. It becomes trivial if any of the 
indices is zero. Applying (<9 2 + d 3 ) ■ k 2 and d 3 ■ (k 3 — k 2 ) to the integrand of (|3.36| ), 
we obtain 



d - n 2 - n 8 - 2n 5 + n 2 2 + + n 8 8 + (7~ 



d — n 3 — n 8 — 2^7 + n 3 3 + l + n 8 8 + (5 



J,, 



0. 



J A = 0. 



(3.54) 
(3.55) 
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When n 3 < and n% 7^ 1, we can raise n 3 by ( [3.48 ); when n 3 < and n 7 7^ 1, we 
can raise n 3 by ( [3.41]) ; when n 3 < and n 8 7^ 1, we can raise n 3 by ( |3.49|) ; when 
n 3 < and n 2 = n 7 = n 8 = 1, we can raise n 3 or n 2 by ( |3.50| ). When n 2 < 
and n 3 7^ 1, we can raise n 2 by Q3.47 ); when n 2 < and n 5 7^ 1, we can raise n 2 
by ( p.41| ); when n 2 < and n 8 7^ 1, we can raise n 2 by ( p.40| ); when n 2 < and 
n 3 — n 5 — n 8 = 1) we can raise n 2 or n 3 by ( 3-50| ). When n 8 < and n 7 7^ 1, we 
can raise n 8 by ( |3.47p ; when n 8 < and n 5 7^ 1, we can raise n 8 by ( 3-48j ); when 
n 8 < and n 3 7^ 1, n 5 = n 7 = 1, we can raise n 8 or n 7 by ( |3.40| ); when n 8 < 
and n 2 7^ 1, n 3 = n 5 = n 7 = 1, we can raise n 8 or n 5 by (3-49|) ; when n 8 < and 
^2 = % = W5 = ?i7 = 1, we can raise n 8 or n 2 by (|3.54|) . When n\ < and n 3 7^ 1, 



we can raise ni by (|3.50| ); when n\ < and 714 7^ 1, we can raise ni using the relation 
symmetric to ( 3-39|) ; when n\ < and n 2 7^ 1, we can raise n\ by ( P-48|) ; when ri\ < 
and n 7 7^ 1, we can raise n\ by ( [3.41D ; when ni < and n 8 7^ 1, we can raise n\ 



by ( |3.49| ); when ri\ < and n 2 = n 3 = n 4 = n 7 = n s = 1, we can raise ni, n 2 or n 3 
by ( |3.42p . When n 5 < and n 8 7^ 1, we can raise n 5 by ( |3.55| ); when n 5 < and 
n 7 7^ 1, we can raise 715 by ( |3.47| ); when ng < and n 2 7^ 1, we can raise or n 8 
by ( 3.49| ); when rig < and n 2 = n 7 = n 8 = 1, we can raise n.5 or n 2 by (|3.48|) . When 



n 7 < and n 8 7^ 1, we can raise n 7 by ( 3.55 ); when n 7 < and 7^ 1, we can raise 
n 7 by ( 3-48|) ; when n 7 < and n 3 ^ 1, n 5 = n 8 = 1, we can raise n 7 or n 8 by (|3.40| ); 
when n 7 < and ni = n 3 = ns = 1, we can raise n 7 or n 3 by (3-47 ). When all the 
indices are positive, we can kill one of the lines 1, 5, 7 by ( 3-55Q . 

Finally we consider J e (ni, n 2 , n 3 , n 5 , n 6 , n 7 , n 8 ) = ld{ni,n 2 ,n 3 ,0,n 5 ,n 6 ,n 7 ,n 8 ) 
(figure |l2|e). This integral vanishes when indices of the following groups of lines are 
non-positive: 6, 13, 18, 17, 58, 78, 57, 325, 237. It becomes trivial if any of the 
indices is zero. Applying (d\ + <9 3 ) • (k 3 — k 2 ) to the integrand of (|3.36| ), we obtain 



d — ni — n 8 — 2n 7 + riil + 3 + n 8 8 + (5~ 



J, ■ = . 



(3.56) 



When n\ < 0, it can be raised using the relation symmetric to (3-51|) . When n 2 < 0, 
it can be raised by ( |3.51| ). When n 3 < and n\ 7^ 1, we can raise n 3 by (3-44 ); when 
n 3 < and n 6 7^ 1, we can raise n 3 using the relation symmetric to ( |3.39| ); when 
n 3 < and n 2 7^ 1, we can raise n 3 using the relation symmetric to ( |3.44j ); when 
n 3 < and n 7 7^ 1, we can raise n 3 by ( |3.39| ); when n 3 < and n 8 7^ 1, we can raise 
n 3 by (| 



; when n 3 < and ni = n 2 = n6 = n 7 



n 8 



1, we can raise n 3 , ni or 



n 2 by (3 .42Q . When n 8 < and n 7 7^ 1 , we can raise n 8 using the relation symmetric 
to (3-48|) ; when n 8 < and n 5 7^ 1, we can raise n 8 using the relation symmetric 
to (3-44 ); when n 8 < and ni 7^ 1, n 5 = n 7 = 1, we can raise n 8 or n 7 using the 
relation symmetric to ( 3.49 ); when n 8 < and n\ — n*, — n 7 — 1, we can raise n 8 
or ni by (|3.56|) . When ra 5 < and n 8 7^ 1, we can raise n 5 by ( |3.56| ); when n 5 < 
and n 7 7^ 1, we can raise n 5 using the relation symmetric to ( 3-50| ); when n 5 < and 
ni 7^ 1, we can raise n 8 or n 7 using the relation symmetric to ( 3-49| ); when n 5 < 
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and rii — n 7 — n s — 1; we can raise ns, n 7 or n\ by ([3.41 ). When n 7 < and 
n$ ^ 1 , we can raise n 7 by ( |3.56|) ; when n 7 < and 715 7^ 1 , we can raise n 7 using the 
relation symmetric to (|3.44|) ; when n 7 < and ni 7^ 1, we can raise n 7 or n 8 using 
the relation symmetric to ( |3.49|) ; when n 7 < and ni = n 5 = n 8 = 1, we can raise 
n 7 or m using the relation symmetric to ( P-48 ). When all the indices are positive, 
we can kill one of the lines 3, 5, 7 by ( |3.56| ). 



3.5 Diagrams with four heavy-quark lines 

Let's define (figure |IT|e) 

N?i 3 N2g 3 d d kid d k 2 d% 



n"i n" 2 n n 3 n n 4 n™5 n™6 n n ? n™8 

-^1 -^2 U Z u <± u h u % u l u % 

= _ i7r 3d/2 ( _ 2a;) 3d +2 (n 1 3+n 2 3)-2Ef =5 ™. ^ ^ n4? ^ ^ n?? ng - ^ ^ 

(& 3 + p) ■ v (h+k 3 +p)-v (k\ +p) ■ v 

L>\ — , LJ2 — , L>3 — , 

D 4 = ^±£tl t D5 = -k h " A, = 

D 7 = -kl^ D s = —{hi — k 2 ) 2 , 



N 13 = 2h ■ k 3 , N 23 = 2k 2 ■ k 3 . 

(3.57) 

This integral vanishes when the indices of the following groups of lines are non- 
positive: 5, 12, 67, 68, 78, 47, 48, 234, 326, 238. The heavy-quark denominators are 
linearly dependent: D\ — D 2 + D 3 = 1, and therefore 

[l-l~ + 2~-3~]ie = 0. (3.58) 

The cases n\ < 0, n 2 < 0, n 3 < reduce to I c , Id- 

Ie{-\n 1 \,n 2 ,n 3 ,n4,n 5 ,n 6 ,n 7 ,n 8 ;n 13 ,n 23 ) = (1 + 1~ - 3~) |ni1 x 

x (6- -4- + 8~) ni3 (5- - 7- + 8- + 0+) n23 / c (n 2 , n 4) n 3 , 0, n 7 , n 5 , n 8 , n 6 ) , 

J e (ni, -|n 2 |,n3,n4,n 5 ,n 6 ,n 7 ,n 8 ;ni3,n23) = 

= (I" + 3- - l)l n2 l(6- - 4 - 8-) ni3 (0 + ) n23 / c (n 1 ,n4,n3,n 5 ,n 7 ,0,n 8 ,n 6 ) , 
h{ni,n 2 , -\n 3 \,n A ,n 5 ,n 6 ,n 7 ,n 8 ;n 13 ,n 23 ) = 

= (1 - I" + 3-) |nsl (4- - 6- + 8- + 0+) ni3 (O+r^K , n 4 , n 2 , n 5 , n 7 , 0, n 6 , n 8 ) . 



If ni,2,3 are all positive, we can lower them by ( 3.58| ) until one of them vanish 
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Let's define (figure |TT]f) 

N? 3 13 N!% 3 d d k 1 d d k 2 d d k 3 

D ni D n 2D n 3D n iD n 5D n 6D nr D n s 



8 

8 



= _, 7r 3d/ 2( _ 2LL , ) 3d +2 (n 1 3 + n 2 3)-2E ! . 5 "*I f ( ni ,n 2 , n 3 , 714, n 5 , n 6 , n 7 , n 8 ; n 13 , n 23 ) , 
(fci+p)-u (fci + fc 3 + p) ■ v (k 2 + h+p)-v 

D A =^2tl t D 5 = -k h " D 6 = -kl, " 

D 7 = -k 2 ? D 8 = -(h-k 2 ) 2 . 

(3.59) 

This integral is mirror-symmetric with respect to 1 <-> 4, 2 <-> 3, 6 <-> 7. It vanishes 
when indices of the following groups of lines are non-positive: 5, 23, 67, 68, 78, 216, 
128, 347, 438. The heavy-quark denominators are linearly dependent: D\ — D 2 + 
D 3 — D4 = 0, and therefore 

[l - - 2" + 3" - 4 - ] I f = 0. (3.60) 

The cases m < 0, n 2 < reduce to J c , i^: 

7/(-|ni|, n 2 , n 3 , n 4 , n 5 , n 6 , ra 7 , n 8 ; n 13 , n 23 ) = (1~ + 2 - 3~) |ni x 

x (6- - 4- + 8-) ni3 (5- - 7- + 8- + 0+) n23 / d (n 2 , n 4) n 3 , 0, n 7 , n 5 , n 8 , n 6 ) , 

i/(ni, -|n 2 |, n 3 , n 4 , n 5 , n 6 , n 7 , n 8 ; n 13 , n 23 ) = (1~ + 2 - 3~) |n2 ' x 

x (4- - 6- + 8- + 0+) ni3 (7- + 8- - 5-) n23 / c (m, n 3 , n 4 , n 6 , 0, n 8 , n 5 , n 7 ) 

(the cases n 4 < 0, n 3 < are symmetric). If «i, 2,3,4 are all positive, we can use ( |3.6C| ) 
to raise, say, n\ and kill one of the lines 2, 3, 4. 
Let's define (figure 

K 3 13 N^ 3 d d k 1 d d k 2 d d k 3 



D^D^Df D2 4 Df Dq 6 DJ} 7 D na 



n n (fci + & 3 + p) • v (fc 2 + &3 + p) • v 

L>\ — , u, l — , L> 3 — , 

Di= (h+P> 1 v t Ds = _ kl " D<j= kl " 

UJ 

D 7 = -k 2 , D 8 = -(ki - k 2 ) 2 . ( 361 ) 

This integral vanishes when the indices of the following groups of lines are non- 
positive: 5, 67, 68, 78, 26, 28, 123, 234, 347, 438. The heavy-quark denominators 
are linearly dependent: D\ — D 3 + D 4 = 1, and therefore 

[l-l~ + 3~-4 - ] J fl = 0. (3.62) 
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The cases n\ < 0, n 3 < 0, n 4 < reduce to lb, Id' 

I g (-\n 1 \,n 2 ,n 3 ,n4,n 5 ,n 6 ,n 7 ,n 8 ;n 13 ,n23) = (1 - 2~ + 3~) |ni1 x 

x (6- - 4- + 8-) ni3 (5- -7- + 8- + + ) n23 I d (n 2 , n 4 , n 3 , 0, n 7 , n 5 , n 8) n 6 ) , 

Ig(ni,n 2 , -\n 3 \,n 4: ,n 5 ,n 6 ,n 7 ,n8;n 13 ,n23) = (1~ + 2" - l)' 713 ' x 
x (4- - 6- + 8- + 0+) ni3 (0+) n23 I d {m, n 4 , n 2 , n 5 , n 7 , 0, n 6 , n 8 ) , 

I g {ni,n 2 ,n 3 , -\n i \,n 5 ,n 6 ,n 7 ,n 8 ;n 13 ,n 23 ) = (1 - IT + 2 - ) 1 ™ 41 x 

x (4- + 6 - 0+) ni3 (4- - 5- + 8-) n23 J 6 (m, n 3) n 2 , n 5 , 0, n 6 , n 8 , n 7 ) . 

If n i,3,4 are an positive, we can lower them by ( |3.62| ) until one of them vanish. 

3.6 Diagrams with five heavy-quark lines 

Let's define (figure \T\\h) 

f N^ 2 N^ 3 N^ 3 d d k 1 d d k 2 d d k 3 = 3d/2( ]3d+ 2( ni2 +n 13 +n 23 )-2Zt 6 ^ x 









x4(m, n 2 , n 3 , n 4 , n 5 , n 6 , n 7 , 


n 8 ,n 12 ,n 13 ,n 23 ) 


D 1 = 


(k 2 + p) ■ v 

•i 


D 2 = 


(ki + k 2 +p) ■ v 


^3 


(ki +p) ■ v 


00 




00 




00 


D 4 = 


{h + k 3 +p) -v 

j 


D 5 = 


(h + p) ■ v 


^6 


- -P 

— ^2 , 




00 




00 






D 7 = 


-P 


D 8 = 


-p 






N 12 = 


2ki ■ k 2 , 


N 13 = 


2fci • k 3 , 


A^ 23 


= 2fc 2 ■ k 3 . 



(3.63) 

This integral is mirror-symmetric with respect to 1 <-> 5, 2 <-> 4, 6 <-> 7. It vanishes 
when the indices of the following groups of lines are non-positive: 6, 7, 8, 12, 45, 
234. There are two linear relations among the heavy-quark denominators: 



[1-1- + 2- - 3-] 4 = 0, [1 - 3- + 4 - 5-] I h = . (3.64) 

The cases 

I h (ni, -\n 2 \,n 3 , -|n 4 |, n 5 , n 6 , n 7 , n 8 ; n 12 , n 13 , n 23 ) = 
= K a (n 1 ,n 5 ,n 3 ,n 6 ,n 7 ,n 8 ; \n 2 \, |n 4 |; n 23 , n 12 , n 13 ) , 

-|n 2 |,n 3 ,n 4 , -|n 5 |, n 6 , n 7 , n 8 ; n 12 , n 13 , n 23 ) = 
= K b (m,n 4: ,n 3 ,n 6 ,n 7 ,n 8 ] \n 2 \, \n 5 \,0,0;n 12 ,n 13 ,n 23 ) , 
I h (ni, -\n 2 \, -|n 3 |,n 4 ,n5,n 6 ,n7,n 8 ;ni 2 ,ni3,n 23 ) = 
= X 6 (ni,n4,n 5 ,n6,n 8 ,n 7 ;0, |n 3 |,0, |n 2 |;n 23 ,ni 3 ,ni 2 ) , 
n 2 , n 3 , n 4 , -|n 5 |, n 6 , n 7 , n 8 ; n 12 , n 13 , n 23 ) = 
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= K c (n 2 ,n 4: ,n 3 ,n 6 ,n 7 ,ns] \ni\, \n 5 \,0;n 12 ,n 13 ,n 2 3) , 
h(-\nx\,n 2 , -\n 3 \,n A ,n 5 ,n 6 ,n 7 ,n 8 ,n 12 ,n 13 ,n 23 ) = 
= K e (n 2 ,n 5 ,n^n 7 ,n 8 ,n 6 ; \nx\,0, \n 3 \;n 12 ,n 23 ,n 13 ) , 

as well as the symmetric cases < 0, rii < 0; n 4 < 0, n 3 < 0; n 5 < 0, n 3 < 0, 
reduce to I c , I d , as will be discussed later in this section. When n 3 < (figure |T2|/), 
we can use 

[1" - 2 + 4 - 5~] 4 = (3.65) 

to raise, say, n\ and kill one of the lines 2, 4, 5. When ni i2 are both positive, we can 
use ( |3.64|) to kill one of the lines 1, 2, 3; when n 4j5 are both positive, we can kill one 



of the lines 3, 4, 5. 

Let's define (figure |i~l"lz) 



/ 



TVlT jV"3 13 jV 2 " 3 23 d d ki d d k 2 d d k 3 _ _■ 3d/2(n \3d+2(n 12 +n 13 +n 23 )-2£f =6 ni 

r-jm n n 2 n n 3 n™* n n s n™6 n«7 n n « ^ ' 

u \ u 2 U Z u <± u % u § u 1 u % 



xli(n 1} n 2} n 3} n 4 , n 5 , n 6 , n 7 , n 8 ; n 12 , n i3 , n 23 ) 



(ki+p)-v (fci + fc 2 + p) • v 

1^1 — , L>2 — , 

UJ UJ 

(fci + k 2 + k 3 + p) ■ v (h + k 3 + p) ■ v (h+p) -v 

U Z — , -^4 — , -^5 — ) 

UJ UJ UJ 

D 6 = -k\ , D 7 = -k\ , D 8 = -k 2 3 . 

(3.66) 

This integral vanishes when the indices of the following groups of lines are non- 
positive: 6, 7, 8, 34, 123. There are two linear relations among the heavy-quark 
denominators: 

[l-l- + 3--4-]ii = 0, [l-l- + 2--5-]/i = 0. (3.67) 



The cases 



Ji(m, -\n 2 \, -|n 3 |,n4,n 5 ,ra 6 ,n 7 ,n 8 ; ^12,^13,^23) = 
= 2f 6 (m, 714,715,716,718,717; |n 2 |,0, \n 3 1 , 0; raw, n 23 , 7li 3 ) 

Ii(-\ni\,n 2 , -|72 3 |,n4,n5,n 6 ,n 7 ,n 8 ;ni2,ni3,n2 3 ) = 
= ^ c (n2,n4,n 5 ,n 6 ,n 8 ,ra 7 ; |ni|,0, \n 3 \;n 12 ,n 23 ,n 13 ) , 

h(-\ni\, -\n 2 \ } n 3} n 4} n 5 ,n 6} n 7} n 8 ;ni 2 ,ni 3} n 23 ) = 
= K d (n 3 ,n 5 ,n^n 7 ,n 6 ,n 8 ; |ni|,0,0, \n 2 \;n 23 ,n 12 ,n 13 ) 
4(-K|,n 2 ,7i 3 , -|n 4 |, n 5 ,n 6 ,n 7 ,n 8 ; n 12 ,n 13 ,n 23 ) = 
= K d (n 3 ,n 5 ,n 2 ,n 7 ,n 8 ,n e ;0, |n, x |, 0, |n 4 |; n 12 , n 23 , n 13 ) 
Ii(ni,n 2 ,n 3 , -|n 4 |, -|n 5 |, ra 6 , n 7 , n 8 ; rii 2 , rii 3 , n 23 ) = 
= K d (n 3 ,ni,n 2 ,n 6 ,n$,n 7 ;0, |n 5 |, |n 4 |, 0; ni 2 , n i3 , n 23 ) 
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Ii{n 1 ,n 2 , -\n 3 \,rii, -\n 5 \, n 6 , n 7 , n 8 ; n 12 , n 13 , n 23 ) = 
= K e (n 4 ,n 1 ,n 2 ,n 6 ,n 7 ,n 8 ;0, \n 3 \, \n 5 \; n 23 , n 13 , n 12 ) , 
h(ni, -\n 2 \,n 3 , -\n 4 \,n 5 ,n e ,n 7 ,n 8 ;n 12 ,n 13 ,n 23 ) = 
= K f (n 1 ,n 5 ,n 3 ,n 6 ,n 7 ,n 8 ; \n 2 \,0, \n 4 \;n 13 ,n 23 ,n 12 ) , 
Ii{-\nx\,n 2 ,n 3 ,n 4 , -\n 5 \,n 6 ,n 7 ,n 8 ; n 12 ,n 13 ,n 23 ) = 
= K g (n 2 ,n 4 ,n 3 ,nQ,n 8 ,n 7 ; |n 5 |, |m|, 0; n 12 , n 23 , n 13 ) 

reduce to I c , Id, as will be discussed later in this section. When n\ < 0, we can use 

[2" - 3 + 4 - 5"] h = (3.68) 

to raise, say, n 2 and kill one of the lines 3, 4, 5. When n 2 $ are both positive, we can 
use ( [3.671) to kill one of the lines 1, 2, 5; when n 3)4 are both positive, we can kill one 
of the lines 1, 3, 4. 

Finally, let's define (figure [TTI7) 

iVT 2 12 iVTo 13 iV 2 n 3 23 d d k 1 d d k 2 d d k 3 



x 



/ N 12 2N l^ N 2i 3d k ld k 2 d k 3 = 3d/2( _ 2 . ) 3d+2(n 12 +n 13 +m 3 )-2Z 8 i=6 n i 

J D" 1 Z>2 2 D 3 3 -D4 4 -D5 5 -Dg 6 D 7 T Z)g 8 V ' 

xlj{ni, n 2 , n 3 , n 4 , n 5 , n 6 , n 7 , n 8 ; n 12 , n 13 , n 23 ) , 
(&i + p) • v (kx + k 2 + p) ■ v 

L>\ — , L) 2 — , 

UJ UJ 

{ki + k 2 + k 3 + p) ■ v (k 2 + k 3 + p) ■ v (k 3 + p) ■ v 
-^3 — j u <± — , -^5 — , 

UJ UJ UJ 

D 6 = -k\ , D 7 = -kl , D 8 = -kj . 

(3.69) 

This integral is mirror-symmetric with respect tol^->5,2^->4,6^8. It vanishes 
when the indices of the following groups of lines are non-positive: 6, 7, 8, 123, 234, 
345. There are two linear relations among the heavy-quark denominators: 

[1 - 1" + 3" - 4-] lj = 0, [1-2- + 3- - 5~] lj=0. (3.70) 

The cases 

Ij(ni, -\n 2 \, -\n 3 \,n 4 ,n 5 ,n 6 ,n 7 ,n 8 ;ni 2 ,n 13 ,n 23 ) = 
= K b (ni, n 4) n 5 ,n6,n7,n 8 ;0,0, \n 3 \, |n 2 |; n 13 , n 23 , n 12 ) , 
Ij(-|ni|, -\n 2 \,n 3 ,n4,n 5 ,ne,n 7 ,ri8;ni 2 ,n 13 ,n 23 ) = 
= K d (n 3 ,n 5 ,n4,n 8 ,n 6 ,n 7 ] \ni\,0, \n 2 \,0;n 23 ,n 13 ,n 12 ) , 
Ij(-\ni\,n 2 , -|n 3 |,n4,n5,n 6 ,n 7 ,n 8 ; n 12 ,n 13 ,n 23 ) = 
= K e (n 2 ,n 5 ,n 4 ,n 8 ,n 7 ,n 6 ; \ni\, \n 3 \, 0; n 12 , n 13 , n 23 ) , 
Ij(ni, -\n 2 \, n 3 , -\n 4 \,n 5 , n 6 , n 7 , n 8 ; n 12 , n 13 , n 23 ) = 
= Kf(m,n 5 ,n 3 ,ne,n 8 ,n 7 ;0, \n 2 \, \n 4 \; n 12 , n 23 , n 13 ) , 
Ij{-\ni\,n 2 ,n 3 ,n 4 , -\n 5 \, n 6 , n 7 , n 8 ; n 12 , n 13 , n 23 ) = 
= -^ 9 (^2,«4,«3,^6,^8,^7;0, |ni|, |n 5 |;ni 2 ,n 2 3,ni 3 ) , 
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a b c 




as well as the symmetric cases n 4 < 0, n 3 < 0; 715 < 0, n 4 < 0; n§ < 0, n 3 < 0, 
reduce to I c , Id, as will be discussed later in this section. When n 3 < 0, we can use 

[1" - 2 + 4 - 5~] Ij = (3.71) 

to raise, say, n\ and kill one of the lines 2, 4, 5. When n 4i4 are both positive, we can 
use ( |3.70| ) to kill one of the lines 1, 3, 4; when n 2j 5 are both positive, we can kill one 
of the lines 2, 3, 5. 

The reduced forms of the integrals Ih,i,j are (figure [TB|) 

K a (n 1 ,n2,n 3 ,n4,n 5 ,n 8 ;n 9 ,n 10 ;n45,n A s,n 5 s) = (IT + 3~ - l) n<J (2~ + 3 - l) ni ° x 
x (0+r 45 (6- - 4 - 8-) n «(7- - 5 - 8-)" S8 J c (n 1 ,n 2 ,n3,n4,n 5 ,0,0,n 8 ) , 

K b (m, n 2 , n 3 , n 4 , n 7 , n 8 ; n 9 , n 10 , n u , n 12 ; n 48 , n 78 , n 47 ) = (1~ + 3~ - l)™ 9 x 
x (2" - 3" + l) ni0 (l - + 2 - l) nil (l" + 2 - 3~) ni2 (6~ - 4 - 8~) n48 x 
x (7- - 5- + 8~)™ 78 (4- -6- + 8- + + )" 47 J c (ni, n 2 , ra 3 , n 4) 0, 0, n 7 , n 8 ) , 

-^ c (ni,n2,n3,n 6 ,n7,n 8 ;n9,nio,nii;n 68 ,n 78 ,n 6 7) = (1~ - 3" + l)™ 9 x 

x (2" - 3" + l) ni0 (l~ + 2- - 3~) nil (6~ - 4 + 8~) n68 (7~ - 5 + 8~) n78 x 
x (4- + 5 - 6 - 7- + 0+)™ 67 7 c (m, n 2 , n 3 , 0, 0, n 6) n r , n 8 ) , 

K d (n 1 ,n 2 , n 3 , n 5 , n 6 , n 7 ; n 9 , ni , n n , n 12 ; n 57 , n 5e , n 67 ) = (l - - 3 + l)™ 9 x 
x (3" - 2 + l) ni0 (l _ -2- + + 2 - 3~)" 12 (5~ + 7" - 8~)™ 57 x 

x (5" - 7" + 8" + + )" 56 (6~ + 7" - 4 - 5 - + )" 67 x 
x J c (m, n 2 , n 3 , 0, n 5 , n 6 , n 7 , 0) , 
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K e (ni,n 2 ,n 3 ,n 5 ,n 6 ,ns]n 9 ,n 1 o,n 11 -,n 68 ,n 5 s,n 56 ) = (1 +2 - 3 )™ 9 x 

x (1~ + 2 - 1)" 10 (3" - 2 + l) nil (6" - 4 + 8~) n68 (7~ - 5 - 8~) n58 x 
x (5" -7- + 8- + + )" 56 / d (ni, n 2 , n 3 , 0, n 5 , n 6 , 0, ra 8 ) , 

Kf(m,n 2 , n 3 , n 4 , n 5 , n 8 ; n 9 , ni , %; n 48 , n 58 , n 45 ) = (l - + 2" - l)™ 9 x 

x (3" - 2 + l) ni °(3~ - 1" + l) nil (4~ - 6 + 8~) n48 (5- - 7~ + 8~) n58 x 
x (0 + )™ 45 / d (n 1 , n 2 , n 3 , n 4 , n 5 , 0, 0, n 8 ) , 

K g (n 1 ,n 2 , n 3 , n 6 , n 7 , n 8 ; n 9 , ni , n u ; n 68 , n 78 , n 67 ) = (1 _ + 2" - 3 _ ) n9 x 

x (3" - 2 + 1)" 10 (3- - 1" + l) nil (6~ - 4 + 8-) rt68 (7- - 5" + 8")™ 78 x 
x (4- + 5- - 6- - 7- + 0+) n67 / d (n 1 , n 2 , n 3 , 0, 0, n 6 , n 7 , n 8 ) . 



4. Implementation and testing 

The package Grinder consists of a set of mutually recursive procedures for Feynman 
integrals of various topologies, which reduce a given Feynman integral to simpler ones, 
until boundary-case integrals with known values are reached. The package is written 



in REDUCE |23], |24|. Remembering results of previous function calls may make the 
program run much faster, if there is enough memory (unfortunately, REDUCE uses 
linear look-up in remember-tables). 



I also re-implemented it in Axiom [25|. All expressions involved are linear com- 



binations of basis integrals with coefficients which are rational functions of d. It is 
convenient to use Axiom domain Vector Fraction UnivariatePolynomial, which 
has all the necessary operations. This makes intermediate expressions shorter than in 
the case when multivariate rational functions are used for entire expressions, because, 
typically, not all basis integrals are accompanied by every possible denominator. The 
amount of GCD calculations is thus reduced. This improvement can be, in princi- 
ple, back-propagated to the REDUCE implementation by using matrices. However, 
working with matrices in REDUCE is awkward, because there are no local matrix 
variables, and no easy way for a function to return a matrix. On the other hand, 
complete diagram calculations, including tensor and 7-matrix algebra, can be done 
in REDUCE. 

The main method of testing was checking various recurrence relations (including 
those which were not directly used for construction of the algorithm) in nested loops 
over rij. For each of two-loop and generalized two- loop integrals, which depend on 5 
indices, a typical number of checks was about 20000; each test set of this size runs 
for a few hours. Integrals with 8 or 9 indices are more difficult to check. Some test 
sets required a few days of CPU time. A tool showing how many times each linear 
code segment has been executed would be invaluable for setting up test cases which 
check all branches at least once. Unfortunately, such a tool is not available in either 
programming system, and I had to emulate it by hand. 
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